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Abstract 

Under certain circumstances, some of which are made explicit here, 
one can deduce bounds on the full sum of a perturbation series of 
a physical quantity by using a variational Borel map on the partial 
series. The method is illustrated by applying it to various examples, 
physical and mathematical. 


1 Introduction 

As there are hardly any problems in physically relevant quantum field theo¬ 
ries that are exactly solvable, one must resort to approximate calculations. 
The only analytical tool that in principle allows for such approximate com¬ 
putations to be systematically improved, and for the errors to be estimated, 
is perturbation theory. Unfortunately in many cases the perturbation pa¬ 
rameter is not small so that the series, of which in practice only a few terms 
are known, diverges or gives a poor representation of the physical quantity. 
Consequently, several techniques have been used to estimate the full sum 
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of an infinite series from the small nnmber of given terms, see for example 
i |,i and references therein. 

A variational method which was introdnced in Ref.|^ has the novelty 
of apparently giving a monotonically converging seqnence of approximations 
and thns possible bonnds on the snm of an inhnite series. That method has 
been applied to the physically important case of the free energy of gange the¬ 
ories at high temperatnre in |^, ^ and the pnrpose of this paper is to provide 
a detailed technical jnstihcation for the methodology nsed there. In essence 
this paper is a condensed and slightly modihed version of the nnpnblished 
report 0 . The exposition here is not meant to satisfy a mathematical physi¬ 
cist bnt is rather snggestive and is illnstrated with several examples to make 
some of the technical assnmptions plansible. It is hoped that some readers 
will hnd the method snfhciently interesting to re-examine it more rigoronsly 
in the fntnre. 

In the following section, I describe in a detailed bnt abstract manner how 
the techniqne of @] can be nsed to obtain lower or npper bonnds for some 
series, and also in many cases to obtain accnrate estimates of the exact snm. 
Concrete examples are given in Sec. (3). A conclnsion and snmmary of open 
problems appears in Sect.(4). Finally, the Appendix contains an approximate 
bnt important analysis of the relevant eqnations which explains some of the 
empirically observed trends. 


2 The Method 


Consider the expansion 


N 


Sivw = i:/„A", 


( 1 ) 


71=0 


with A the pertnrbation parameter. It will be assnmed in this paper that 
the relevant expansions are inhnitely long, bnt that only a hnite nnmber of 
terms are known. Since pertnrbation expansions in qnantnm held theories 
have the generic behavionr /„ ~ n\ for n large []^, it is natnral to consider 
the Borel transform 



( 2 ) 


The series (0) can then be recovered from the Borel integral. 



( 3 ) 
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If the exact function B{z) = B^{z) is known, Eq.(^) can be taken as defining 
the exact sum, S = Soo, of the full perturbation series if the Borel integral is 
well-dehned. The poor convergence of can be attribnted to the expansion 
of B{z) in a power series as in (|j) and the snbseqnent nse of that series 
beyond its radins of convergence in (|^). 

Snppose initially that B{z) has only one singnlarity in the complex 2 ;- 
plane (the Borel plane), at 2 ; = —1/p, with p real and positive. Then the 
radius of convergence of the Borel series (^) is 1/p. Therefore in order to 
constrnct an approximation to the exact snm, one must extend the domain 
of convergence of the partial series in Eq.(^. The method of Loeffel, Le- 
Gnillon and Zinn-Justin 0,i is to use a conformal map 


w{z) 


\/l + zp - 1 
a/I + zp + 1 ’ 


( 4 ) 


which maps the 2 ;-plane into a nnit circle in the w-plane, with the singularity 
at z = — 1/p mapped to tc = — 1. The inverse of this map is 


z{w) 


Aw 1 
p (1 — wy 


Making the change of variables in (|^) one obtains. 


( 5 ) 


^ 1 /*! fj y 

S^{\) = - dw — B^izH) , (6) 

A Jo dw 

where B]^{z{w)) is understood as a power series in w obtained by an expan¬ 
sion of (^. In terms of the variable tc, the series Bn{z{w)) converges for 
|w| < 1, so that the potential problem in (^) is only at the npper limit of 
integration. Now, the difference between Bn{z{w)) and Bi\f+i{z{w)) begins 
at order w^^^. If one knows the coefficients /„ only up to n = N, then it is 
consistent to keep only terms np to order w^ in Bi^{z{w)). Performing this 
trnncation in (|^) and riverting back to the variable, one obtains 


^v(A) 


1 

A 


n=N 


E 

n=0 


fn 


n\ 


/4y {2n + k-l)\ 

U/ ho k\{2n-l)\ 



dz . (7) 


SnW is a non-trivial resnmmation of the original series In Ref.|||, 

0, Eq.(^ has been used to resum long pertnrbation expansions for critical 
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exponents, with z = —\/p the location of the instanton singularity of held 
theory. Note that after a scaling, m may be written as 


5iv(A) 


n=N 


E 

n=0 


fn 


Ay y." (2n + fc- 1)! 
U/ to k\{2n-l)\ 



dz e-" 


( 8 ) 


Since w{z) is a bounded and slowly varying function, this shows that the 
resummed expression Sn{^) is a much slower varying function of the coupling 
than the original series Sat (A). 

The resummation (J^) has also been applied to quantum chromodynamics 
(QCD), with z = —1/p the location of the hrst ultraviolet renormalon pole 
Bill. Actually, QCD is Borel-nonsummable [0, which means that B{z) 
also has poles for z > 0, rendering the Borel integral ambiguous. In such a 
case, one can still use (|]) to dehne the sum of the perturbation series once 
the integral is made dehnite through some prescription such as the principal 
value. 

In all the applications of in [|, H, |^, p is a known and hxed 
constant which determines the location of the singularity of B{z) closest to 
the origin. In Ref. |^, the expression (|^) was used as the starting point for a 
new technique which can be used even in cases when the singularity structure 
of B{z) is unknown. Indeed, as seen above, the conventional transformation 
of a series 


N 


Sn{\) = E 


(9) 


n=0 


by the Borel-conformal map method results in the reorganised expression 


-| n=N 

Sn(\p) = y E 

n=0 


fn {2n + k-l)\ 

k=0 


, dz w(zY^+^\ 
\P J S k\{2n — 1)\ Jo 

( 10 ) 

with p a fixed constant. However, notice that p does not hgure in (|]) but 
enters (|^ only through the conformal map (^. Thus instead of treating p 
as some hxed constant as in Ref.lSl, 


10 , 11 


one may consider p > 0 a free 
parameter which dehnes the conformal map and Eq. (^) then represents 
a continuous family of resummations, one for each value of p. Thus, from 
now on, p will not refer to the location of some singularity in B{z). In fact, 
except for the regularity of B [z) at the origin, no other knowledge about the 
singularity structure of B{z) is required for the method elaborated below. 
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Of course having liberated ourselves from the usual interpretation of p 
in ( 0 ), we also lose dehniteness in our resummation. Therefore some new 
condition must be imposed to £x the value of p and hence of the expression 
( 0 ). For a start, for each iV, choose p > 0 to be the location of an extremum 
of Sn{^,p) 0 . Since ([T^) depends on A, the value of p in principle will also 
depend on A. However the procedure then becomes too unwieldy, and to 
simplify it p is determined at some reference value A = Aq, say at the mid¬ 
point of the range of interest: 


g5'jv(Ao,p) 

dp 

As mentioned after (H), S{X,p) is a relatively slowly varying function of A, 
hence determining p at some hxed A = Aq, and then using the same p in (|IoD 
for various A is sufficient for practical purposes. Indeed, in many applications, 
one requires the sum of the series at one particular value of the coupling, or 
in a very narrow range, so the simplihcation made in (PH]), is both useful and 
sufficient. 

In some cases, (^) will not have any solutions. For example, if all the /„ 
are of the same sign, Sn{X,p) will be a monotonic function of p. Thus for 
the procedure to work, at least some of the /„ must be of a different sign. 
This will be assumed to be the case from now on. Suppose next that /i 7 ^ 0. 
Then for p large, and since the p dependence of w{z) is mild. 



^iv(p)~/o + -. (12) 

p 

Therefore for /i < 0, will first decrease as 1/p increases and then increase 
when the next term fm/p^ > 0 dominates the sum. If the last non-zero 
term /n/p^ is positive then one expects Sn{p) to have a global minimun at 
some p > 0. However if Jn/p^ < 0 then Sisi{p) can become very small as 
p —O’*" and so the minimum would likely be only a local extremum. From 
this heuristic argument one concludes that if the hrst non-zero coefficient 
fn, {n > 0) is negative, Eq.([Tl|) will have global minima solutions for some 
N. 

For conciseness, unless otherwise stated, I will discuss from now on only 
the case when a global minimum exists, as the arguments for the global 
maximum case are then obvious. For a series S'jv(A,p), define p{N) to be 
the position of the global minima for the values of N when they exist, and for 
other values of N let p{N) denote the location of the local minima. Also, let 
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Sn = Sn{X,p{N)). (13) 

Obviously the definition of p{N) has been chosen because it is useful. If 
p{N) is a location of a global minimum, then certainly is a lower bound on 
Sn{X,p). However this does not imply that Sjy is a lower bound on the exact 
sum S of the pertutbation series because for finite N, S might he outside of 
the space of resummations labelled by p. For each N, let p*{N) be the value 
of p that is optimal, that is, it is the value which when used in (| 10 |) gives the 
best estimate of S. Define, = Sn{X,p*{N)). Then for a global minima 
one has 

Sn < (14) 

Presumably converges to S' as iV —> cx). Then for those N when global 

minima exist, 

Sn^oo < s. (15) 

(This implicitly assumes that the sub-sequence of global minima is infinite: 
That is, given any positive integer Nq, there is some n > Nq for which /„ is 
positive.) 

Though the inclusion of global minima in the definition of p{N) is quite 
clear, the inclusion of local minima requires some explanation. As mentioned 
earlier, if Sn{X.p) has a global minimum, then Sn+i{X,p) will probably not 
have a global minimum if /n+i is negative because Sn+i{X,p) might become 
very small as p ^ 0"^. However S'Ar_|_i(A,p) will still have a local minimum 
(and hence also a local maximum). The local minimum will occur for mod¬ 
erate values of p close to p{N), the global minimum position of Sn{X,p). 
Therefore one might expect that the local minimum for S]s[+i{X,p) still gives 
a bound on Indeed, suppose that the inequality 

Sn < Sn+1 (16) 

holds for all N up to = oo, then clearly 

Sn<S (17) 

and one concludes that Sn{X,p{N)) is a lower bound on the sum of the 
full perturbation series. Of course proving (|I6|) is equivalent to knowing /„ 
explicitly for all n, which is not the situation in reality. As a practical matter, 
it is sufficient to draw a conclusion after observing the trend ([l 6 |) (or lack 
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of) for the available terms of the series. Arguments given in the Appendix 
indicate that the trend, once started, will continue. 

An interesting question is whether the inequality (pA]) is saturated. In 
most of the examples studied this has been found to be the case. In fact 
the convergence is so rapid that one can conclude with a high degree of 
conhdence, not only that is a. lower bound, but that it is close to the 
exact value. However the example in Sec. (3.5) shows that the series S^, 
though forming a bound on the actual value S, might not converge to S. 
In that example the exact result ^(A) has a hrst derivative dS{X)/dX that 
varies rapidly with A. An explanation of why in such cases the bounds Sn 
do not converge to the exact value is given in the Appendix. However, as 
described below, such situations can also be taken care of by the resummation 
procedure (p!0| - |TT|) . 

Recall that when a local minima occurs for some N, one expects also a 
local maximum. Define pN to be the position of those local maxima and Sn 
the corresponding value of the resummed series. It turns out that Sn also 
satishes an inequality like ([T^ ) (see Sec. (3.5)). However since for the sequence 
S'tv of only local extrema one does not seem to have a statement like (|T^|T^, 
it is not a priori obvious that they form bounds to the exact result. Of course 
if one believes that the p*{N) as dehned above always take moderate values, 
then a statement like (^) might also be made for the sequence Sn- 

Furthermore, given a series which has global minima solutions to (pd^, 
then by ’removing’ the hrst nontrivial coefficient of that series one forms an 
auxiliary series, 5^(A, p), which has global maxima solutions. In this way one 
can form alternative bounds to the sum of the original series complementing 
those obtained from Sn- This is illustrated in Sec.(3.5). 

Using the sequences Sn, Sn, and 5^, one can obtain constraints on S 
and, with some physical input, an estimate of S itself. A question now arises 
for practical problems where exact results are not known. Is it possible to 
determine whether the unknown exact result lies closer to the upper or lower 
bound ? Empirically, it appears that the series Sn converges to the exact 
result ^(A) if d^S{X)/dX^ is very small in magnitude in the entire range of 
interest, 0 < A < Aq. Otherwise the exact result seems to be better approxi¬ 
mated by the series Aat or by the auxiliary series method discussed above. An 
explanation of why the bounds Sn have slowly varying hrst derivatives, and 
why the bounds Sn and those from the auxiliary series have faster varying 
hrst derivatives is given in the Appendix. 

Although the main discussion in this paper will be for the p{N) (or p{N) 
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andp'(iV)) as defined above, in some cases one finds (by inspection) solntions 
Po{N) to (0) which have the property that as iV — oo, po{N) pq, a. 
constant. In snch a case one is tempted to speculate that the hxed point po 
indicates the existence of a singularity at z = —1/p in B{z). This is indeed 
found to be the case in the examples studied, though apparently the converse 
is not necessarily true: the singularities of B{z) need not show up as solutions 
of (0. Furthermore, the convergence of the series S'(o) 7 v is not expected to 
be monotonic since in general the po{N) refer to both maxima and minima. 

Though not manifest at hrst sight, the analysis in the Appendix suggests 
that even the p{N) dehned above Eq.(0) actually will converge to a hxed 
value as N ^ oo. This trend can be observed in the examples studied. 


3 Examples 


In this section a number of examples with different properties are studied 
to illustrate the resummation technique of Eqs.(pT|-|TT|) and to show that the 
various technical assumptions made there are not vacuous. Except for the last 
example, all dehne Borel-summable series, that is, B{z) has no singularities 
at real, positive 2 ;. A Borel-nonsummable example will be considered in 
Sec.(3.3). All of the examples here have global minima as solutions to the 
extremum condition (0) for some N. (From any such series Sn one can 
trivially construct C — Sn, with C a constant, which then gives an example 
of a series with global maxima solutions to (0)). The reader is reminded 
that the existence of global minima does not by itself imply that one has 
obtained a lower bound on the sum of the series. The additional condition 


(0) must be satished for a lower bound. The Euler-Heisenberg series below 
shows how one can end up with an upper bound from global minima ! (See 
also the Appendix). 

It is also re-emphasized that although in these examples the exact singu¬ 
larity structure of B{z) is known, that information will not be used in the 
resummation. The resummation of the partial series 


N 


Sn{X) = 1: fnX 


(18) 


71=0 


will proceed using Eqs.(0-0). The exact B{z) will only be used to compare 
the resummed results with the exact sum of the full series given by the Borel 
integral 




1 

S(X) = - e-^/^Biz)dz. (19) 

A JO 

In addition to the examples below, the interested reader may hnd others 


m 


3.1 A Prototype 

Consider a model defined by the Borel fnnction, 

1 


B(z) = 


1 + z 


Expanding B{z) to A^-th order in z and using the result in 
truncated series 


N 


sn = 


n\ 


( 20 ) 

gives the 

( 21 ) 


n=0 


The divergent nature of this series is displayed in Fig.(la). Using /„ = 
in (pUD and solving Eq.([TT|) at the reference value Aq = 1 gives 
the following solutions (minima): p(2) = 2.65, p(3) = 5.1 and p(4) = 8.4. 
As expected from the arguments of Sec. (2), there is no solution for N = 
1, the solutions for N = 2 and = 4 are global minima while that for 
A^ = 3 is a local minimum. The resummed series is shown in Fig.(lb). The 
convergence of the Aat is monotonic and satishes the condition (p!6|). Hence 
the approximants Sn can be argued to form lower bounds to the exact result. 
That this is indeed the case can be seen by inspection of the exact result in 
Fig.(lb) as obtained from (|^) and (0). Furthermore, it is clear that the 
lower bounds converge to the exact value, and so may be used to estimate 
it. At A = 0.5, the exact value is 0.722657, while the approximants are 
S 2 = 0.704, S 3 = 0.709, S 4 = 0.711. 

For this model there are other solutions, for each N, to the extremum 
condition (^) in addition to the minima. By inspection one picks out the 
sequence of values, Po{3) = 1.6 (local maximum), po(4) = 1-3 (local mini¬ 
mum) and po(5) = 1.15 (local maximum) as plausibly approaching a hxed 
point. Indeed we already know that in this model the exact singularity of 
B{z) is a pole at z = —1/p = —1, so the hxed point probably refers to the 
location of this singularity. It should come as no surprise that if the approx¬ 
imants (1^) are evaluated at the values po{N), the convergence to the exact 
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value will be much faster, and this is indicated in Fig.(lc). At A = 0.5 the 
values of the approximants are S(o )3 = 0.726, S'(o )4 = 0.7219, S(o )5 = 0.7228. 
Notice however that in this case the convergence is not monotonic (and was 
not expected to be). 


3.2 Beta function of a Two-Dimensional Field Theory 


In [131 the exact beta function in a two-dimensional held theory was obtained 


and the result checked explicitly up to three-loop order. After a scaling, that 
beta function is essentially of the form 


^(A) = 


1 + A 


( 22 ) 


In this case a power expansion of S(A) is actually convergent for |A| < 1. 
However close to A = 1, the convergence is very slow and one reqnires a 
large number of terms of the series to obtain accurate results. The truncated 


pertnrbation series corresponding to (^ 2 ]) is 


N 

= ( 23 ) 

0 

Notice that there is no factorial growth of the coefficients in (|23|) . The slow 
convergence of this series for A close to 1 is displayed in Fig. (2a). 

The ntility of the resnmmation procedure in this case is now demon¬ 
strated. Using fn = (—1)*^ in (1T0|) and solving Eq.(|Tl|) at the reference 
value Ao = 1, gives the following solutions (minima): p{ 2 ) = 1 . 2 , p(3) = 
2.9, p(4) = 5. The resummed series is shown in Fig.(2b) together with the 
exact resnlt. The convergence is rapid, monotonic, satishes the condition 
(H), and the approximants Sn form lower bounds to the exact resnlt. With 
only fonr terms, the resnmmed series already shows a dramatically improved 
convergence compared to the original series (^3]) , even for conplings as large 
as A = 0.8: The exact snm at that conpling is 0.556, while the resummed 
values are S 2 = 0.492, S '3 = 0.507, S 4 = 0.512. 

This example illustrates that the resummation procednre (pUj - pd]) is usefnl 
even for a convergent series, especially when one is close to the radius of 
convergence and when only a few terms of the series are available (which is 
often the sitnation in practice). Furthermore this example emphasizes that 
the parameter p has no obvious relation to the singularities of B{z). Indeed 
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the function (P^ ) corresponds to the Borel transform 

B{z) = e-" . (24) 

which is regular everywhere in the Borel plane. 


3.3 A Borel-Nonsummable Model 


For most of the physical quantities calculable from the Standard Model of 
particle physics, the perturbation expansion is not expected to be Borel 
summable. In simple terms, this means that the function B{z) has poles 
on the positive semi-axis of the Borel plane thus rendering the Borel integral 
(|T^) ambiguous. If one choses an ie presciption then the resulting ambigu¬ 
ity is in the imaginary part. Sometimes these imaginary parts are of direct 
physical relevance [|I4|, More generally they indicate that the perturbation 


series does not give the full answer, but must be supplemented with some 

0 . 


non-perturbative contributions |10 

,-l/A 


Since the imaginary parts will be of 


the form the additional real non-perturbative terms are expected to 

take the same form. These expectations have been conhrmed in some lower 
dimensional held theories (see 


1 ^ and references therein). 


One can also construct mathematical models to illustrate the arguments 
of the last paragraph. Suppose, for simplicity, that the only singularity of 
B{z) for z > 0 is a single pole at z = q. Then the sum of the perturbation 
series can be dehned by the principal value prescription. 


1 

Sp,rt = ^V dz B{z ). (25) 

A Jo 

With this dehnition, one focuses only on the real part of the physical quantity. 
Suppose furthermore, again for simplicity, that B{z) is integrable at inhnity. 
Then the Eq.(|25|) may be rewritten as 


1 roo p~qI^ roo 

Spert = Y / B{z) + / dz B{z) . (26) 

A Jo A Jo 

The hrst term on the right-hand-side is hnite and unambiguous. Call it Sexact, 
and denote the second term on the right-hand-side as Snp- Thus we have 

SexactiX) = SpertiX) - ^np(A) . (27) 
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In this way, the exact result Sexact, has been broken into two components, 
Spert which is purely perturbative, and Snp which is purely non-perturbative. 
However while Sexact is well-dehned, both Spert and Snp are only dehned 
through the principal value prescription. Though highly simplihed, this 
model plausibly represents the situation, for example, in Quantum Chro¬ 
modynamics (QCD). 

In order to illustrate explicitly the resummation technique ([IO| - pdD in the 
Borel-nonsummable case, set 


B{z) = 


Then, 


5'peri(A) = jV 



(28) 

—zjX 

z){h-z) ■ 

(29) 


The truncated perturbation series corresponding to 


IS 


N 




(30) 


with 

/„ = ((-l)" + 5-("+')) (31) 

0 

The divergent series is oscillatory, analaogous to that Fig.(la). The solution 
of Eq.(pd]) at the reference value Aq = 1, gives (minima): p(2) = 2.8, p(3) = 
5.4, p(4) = 9, p(5) = 13.5. The resummed partial series is shown in Fig.(3a) 
together with the exact sum of the full perturbation series given by (^9]). 
The convergence is rapid, monotonic, satishes the condition ([T6|), and the 
approximants S'tv do indeed form lower bounds. Also, the convergence is so 
fast that, except for a small interval at intermediate coupling, S^ is not only 
a lower bound, but also a very good approximation to Spert- 

Now, corresponding to the Borel function (P8[), one has from the dehni- 


tions before Eq.(|^), 

^np(A) 

- e-5/A 

6 A 

(32) 

and 

1 

roo p~^l^ _ p~^/^ 




Jo {l + z){5 — z) 

(33) 
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In Fig. (3b), the curves for Sexact and S '5 are plotted for a bigger range of 
A. This shows that except for a small range of couplings, the resummed 
perturbative approximant S 5 lies above and deviates signihcantly from the 
exact result although, as seen above, S '5 does form a converging lower bound 
to the perturbative component Spert of the exact result. The difference is 
of course the non-perturbative piece Snp- In the same hgure there is a plot 
of (S '5 — Snp) which, according to the dehnition (^Tj) should approximate 
Sexact- Indeed the agreement is very good and, amusingly, it improves at 
large coupling: At A = 5, Sexact = 0.0533 and S '5 — Snp = 0.0457, while at 
A = 10 

; Sexact 0.0219 and ^5 - Snp = 0.0193. 

This toy model discussion illustrates concretely the arguments given in 
Ref. [Q for the free-energy density of thermal SU{3) gauge theory. There it 
was found that the exact result, given by lattice data, differed signihcantly 
from the resummed perturbative result, and it was argued that the differ¬ 
ence was caused by the Borel-nonsummability of the theory. Assuming a 
non-perturbative component of the form y the constants A and q 

were determined from the difference between Sexact and S 5 . A more detailed 
discussion of the results in and their extension to full QCD and other 
non-Abelian gauge theories is presented in p. 


3.4 The Euler-Heisenberg Series 


Schwinger’s effective Lagrangian (density) for QED in a uniform mag¬ 


netic held B is. 


L(A) 


^00 (ig 

Jo 



(34) 


where A = e^B'^/nrA, is the renormalized electron coupling, and m the 
electron mass. 

Dehne the dimensionless quantity 


S(A) 



(35) 


which is related in an obvious way to (^). An expansion of (^) is given by 


5(A) = 1600 ^/„A" (36) 

n=l 
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with 


in—1 


fn = 


B- 


'2n+2 


(37) 


2n (2n + 1) (2ij + 2) ’ 

and where B 2 n are the Bernoulli numbers which alternate in sign, and diverge 
factorially for large n. Equation (^) is the Euler-Heisenberg [16| expansion 


which is equivalent to a sum of an inhnite number of Feynman diagrams con¬ 
sisting of one closed fermion loop with an even number of external photons. 
The Euler-Heisenberg series diverges for large values of A and displays oscil¬ 
latory behaviour analogous to that in Fig.(la). Consider now a resummation 
of the divergent series using ( 0 - 0 )- At the reference value A = 10, the 
solutions to ( 1Tl|) are (global minima) p{2) = 0.77, p(3) = 1.4, p( 6 ) = 2.25. 
No solutions were found for N = 3,5, 7. 


The exact result of Schwinger, given by (|35|) , is plotted in Fig. (4), together 
with the approximants S'2, S'4, Sq. The approximants form upper bounds 
although the p{N) are positions of global minima. See again the caution 
prior to Eq. dH) and the Appendix. The convergence of the approximants 
Sm to the exact expression is manifest. For example, at A = 10, where 
the Euler-Heisenberg series is badly divergent, the Schwinger’s exact value is 
5(10) = —8.056 while the estimates are, S 2 = —5.9, ^4 = —6.9, Sq = —7.3. 

The Euler-Heisenberg series (^) is Borel summable |l^, and therefore 
provides an example of a Borel-summable series in a theory (QED) which is 
generally considered Borel-nonsummable. However it should be noted that 
the Euler-Heisenberg series (^) represents a one-fermion-loop result whereas 
the ’’renormalon” singularities [T^ which signal Borel-nonsummability are 
expected to arise when multi-fermion loop corrections to (|3^) are taken into 
account. 


3.5 Critical Exponents in Three Dimensions 

The three-dimensional 0{M) symmetric held theory can be used as an 
effective theory (see [§, || and references therein) to describe the critical 
behaviour of many physical systems near a second-order phase transition. 
For this purpose, the renormalization group functions of this theory have 
been calculated to very high order. A detailed list of references can be found 
in 1^, ^|. For example, the beta function for the polymer case. A/” = 0, is 
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given by 

/3(A) = -A+A2-0.439815A3+0.389923A^-0.447316A®+0.633855A®-1.03493A^ , 

(38) 

where A is a dimensionless coupling (usually denoted as g or g in the literature 
1^). The objective is to hnd the non-trivial infrared hxed point, A*, of the 
theory, which is given by the zero of the beta function, 

= 0, (39) 

= 0^ > 0 . (40) 

The expression ( |38| ) is divergent for large A ~ 1 where a nontrivial zero 
is expected. Using the resummation (0) with Sn denoting approximants 
to the beta function, and choosing the reference point Aq = 1 , one gets as 
solutions to Eg . ([Til) (minima), p(2) = 1.3, p(3) = 3.2, p(4) = 5.6, p(5) = 
8 .6, p{6) = 12.25, p{7) = 16.5. The curves are shown in Fig.(5a). They form 
rapidly converging lower bounds but intercept the A-axis only at the origin, 
thus giving only a trivial zero to the beta function. The situation is similar 
to the sin(7rA) toy model studied earlier. 

Now, since for N = 3, 5, 7 the minima are only local, one expects lo¬ 
cal maxima to exist also. Indeed the local maxima are located at, p{3) = 
0.18, p(5) = 0.21, p(7) = 0.19. The curves for Sn are shown in Fig(5b). 
Though these are local maxima, they obey the inequality Sn+i > Sn and so 
appear to form lower bounds ! Since the lower bounds due to Sn are higher 
than the lower bounds due to S'at (Fig. (5a)), one may argue that those due 
to Sn provide more accurate information. The curves in Fig. (5b) do in fact 
have a non-trivial zero. For N = 7, the zero is near A = 1.425. It is reason¬ 
able to suppose that as both the approximants Sn and S{N) provide lower 
bounds, it is simply a matter of choosing the highest lower bound to estimate 
the sum of the series. 

Given the physical importance of this example, it is useful to perform 
further checks. So let us re-analyse the problem using the auxiliary series 
method discussed in Sect.(2). That is, let Sp^{X,p) = Sn{X,p)/ X. The 
auxiliary series will then obviously have maxima as solutions to ([TT|). The 
solutions at the reference point Aq = 1 arep'(2) = 0.6, p'{3) = 0.875, p'(4) = 
0.33, p'{5) = 0.4, p'{6) = 0.2525. (Actually, p'(5) is a point of inflexion). The 
curves for AS*^ are shown in Fig. (5c). Again, though determined by positions 


/3(A*) 

d(5 
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of maxima, the curves appear to form rapidly converging lower bounds which 
have a non-trivial zero. 

In order to compare with the results in the literature, more precise num¬ 
bers are now quoted. Conjecturing that the curves in Fig. (5c) are indeed 
lower bounds to the exact result, the N = 6 curve, shown magnihed near 
its nontrivial zero in Fig.(5d), gives an upper bound of 1.4193 to the non¬ 
trivial zero of the beta function. Re-optimizing the N = 6 equation (^) at 
Ao = 1.419 does not change the curve for XSj^ signihcantly to modify that 
bound at the level of accuracy quoted. The slope of the beta function at the 
non-trivial zero can also be determined from the N = 6 curve in Fig.(5d). 
It is a; = 0.7955. Since the curves appear to get steeper as N increases, this 
value of a; is a lower bound. These values can now be compared with those 
of Ref. 0: there it was found that A* = 1.413 ± 0.006 and u = 0.812 ± 0.016. 
The agreement with the bounds found here is excellent. (Comparison of 
results obtained by other methods and other authors may be found in |Q). 

It is important to note the difference between the methodology used in 
this paper and that employed in Ref. 0. In 0 the authors used a Borel-Leroy 
transformation with a variable parameter b but they used the conventional 
conformal map with the value of p fixed at the precise location of the 
instanton singularity, p* = 0.166246 (the value for A/” = 0). The parameter 
b, together with some other parameters introduced in |Q] were used to check 
the convergence of the series and to estimate their errors. Here instead 
the usual Borel tranform (|^) is used but the conformal map has a single 
variational parameter p determined according to the condition ([TT|). From 
the numbers quoted above, it is clear that the values p{N) used here are not 
the same as the value p*. Furthermore, in the approach of this paper, the only 
assumptions made about the analyticity structure of the Borel transform are 
(i) B{z) should be regular at the origin, and (ii) in order for the resummed 
perturbation series to faithfully represent the physical quantity, the series 
should be Borel summable. Borel summability of the 
established in 


I theory has been 


From this example it is clear that the novel resummation presented here, 
with the parameter p determined from ([TI|) , can be used to complement the 
analysis done in P]. 
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4 Conclusion 


The method used in Ref. [Q to resum the perturbative free energy of gauge 
theories at high temperature has been detailed. It can be used more generally 
to obtain bounds on the full pertubation expression, S'(A), of a physical 
quantity if in addition to the given partial perturbation series 


N 


Sn{X) = E 


(41) 


n=0 


some additional requirements are met and some assumptions made. 

The first requirement is that not all the /„ be of the same sign. Then 
solutions exist to the extremum condition (|^. If the hrst nontrivial co¬ 
efficient fn{n > 0) is negative, then the solutions to o will be global or 
local minima, depending on the sign of /at. The approximants Sat as dehned 
through then form lower bounds to ^(A) if the inequality (pffiD is satished 
for all N and if the assumptions leading to / 1I3|) are admitted. The inequality 
([T6| ) must hrst be explicitly tested for the available terms of the series, then 
the arguments given in the Appendix suggest that the trend will continue 
for those cases where c{N) > 1 at low orders. Therefore even with partial 
information as in one can apparently deduce lower bounds to the exact 
value 5'(A) if the technical assumptions mentioned above are valid. Under 
similar conditions, additional bounds may be obtained by using the auxiliary 
series method described in Sec. (2). 

Sometimes one hnds minima solutions to (ED but for which the Sn obey 
an inequality opposite to (|^. In such cases, even though the convergence is 
monotonic, it is not a priori obvious that the Sn are actually upper bounds to 
the exact value. For problems where the exact result is unknown, additional 
input from theory or physics is required before a dehnitive statement can 
be made. However in all the examples encountered of this type, the Sn did 
actually bound the exact result. A similar caveat concerns the approximants 
S defined near the end of Sec. (2). 

The observed rapid convergence of the bounds has been explained in the 
Appendix. The obvious question is whether the bounds converge to the exact 
value itself. Empirically, it is found that the bounds formed hj Sn actually 
converge to the exact value if *S'(A) has a slowly varying hrst derivative, 
dS{\)/d\. Otherwise the complementary bounds formed by 5(A) or the 
auxiliary series method give better approximations to the exact result. The 
reason for this phenomena has been given in the Appendix. 
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Although most of the observed features of the resummation method devel¬ 
oped here have been explained in the Appendix, at least in a semi-quantitative 
way, more patterns were detected than could be explained. For example, (i) 
Is it true that in all cases, as N —> oo, c{N) = p{N +1)/p{N) —> 1 (Similarly 
for the p{N)) ? (ii) Is it true that in all cases the approximants Sn and Sn 
form bounds to the full perturbative result ? (hi) Is it true that the bounds 
Sjsf always approximate the exact result S well if I is small through¬ 

out the range, 0 < A < Aq, of interest, and otherwise the alternative bounds 
Sn or those from the auxiliary series give better approximations ? 

Clearly it would be desirable to have a better understanding of the con¬ 
ditions under which some of the technical assumptions made above hold 
rigorously. Until that is achieved, and in particular for those physical cases 
where little is known about the exact series anyway, the method can neverthe¬ 
less be used as another independent estimator of the sum of a perturbation 
series, complementing other techniques in the literature. I point out that 
the resummation method discussed here can be viewed as an example of an 
’’order-dependent mapping” introduced in |T^ by Seznac and Zinn-Justin, 
but applied here to the Borel series rather than the original perturbation 
expansion. In addition, the extremum condition ( pA]) can be thought of as 
an example of Stevensen’s Principal of Minimum Sensitivity . 


Acknowledgements: I am grateful to C. Coriano', S. Pola, J.S. Prakash, 
A. Goldhaber, P. Van Nieuwenhuizen, I. Parwani and U. Parwani for hospi¬ 
tality during the course of this work. 
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Appendix 

A1 

The resummed perturbation series with a variable parameter p is given 
by 


Sn{\p) 


n=N 


E 

n=0 


fn 


uy {2n + k-l)\ 
Uy to k\{2n-l)\ 



dz e-^ (42) 


Write the equation above in the compact form 


Sn{\p) - 2 ^ 

n=0 P 


(43) 


where An{N) is a p and A dependent constant that can be read off from (^21). 
In particular note that the sign of An{N) is the same as the sign of fn- The 
extremum condition (pT|) applied to (^3]) results in 


OSm 

dX 


^ nAnjN) 

hi pW" 


(44) 


where use has been made of the form w{zX) in transforming the p derivative 
of An{N) into a A derivative. At the solution p = p{N) of (pH]), one has 


Sn = 


^ An{N) 

ho p(E" 


(45) 


Since for a given N there is in general more than one solution to the 
extremum equation (^), here p{N + 1) and p{N) will refer to solutions at 
consequtive orders which are both positions of minima or both positions of 
maxima. Now dehne 

p{N + 1) = c{N) p{N) (46) 

where c{N) is some function of N. Then from equation (|4^) and the corre¬ 
sponding one at order -|- 1, one easily deduces 


^d{SN+i-SN), ^ n (An{N + l) \ (iV + l)A^+i(iV + l) 

dX hiP{NY\ Y c(Ar)^+ip(iV)^+i 

(47) 

Now for large N, An{N) ~ An{N + 1), and write this simply as Also 
dehne, AS'jv = — S^. Assume now that p{N) is large and c{N) > 1. 
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Also assume that all the coefficients are generically of the same order, or 
at least do not increase rapidly with n (the factorial growth of /„ has already 
been taken care of by the Borel transform). Then keeping terms needed to 
solve for l/c(iV) at leading order, (^Tj) simplihes to 


A 


dX 




1 

P(.N) 

\c{N) ) 


2 A 2 
p{Ny ■ 


(48) 


Similarly, from eq. (|45|) and the corresponding one at order A^ + 1 one deduces 
at large N and to leading order in l/p{N), 


ASn = 





(49) 


where it is implicit in this equation that A = Aq. Now, if Aq is varied, then 
the leading change in (|49D comes from Ai. Thus using (^), Eq.(|48|) can be 
simplihed to 

—^ (— - 1 

p{N) \c{N) 

where a = {Ai — A^)|ao- Thus, 


2A2 

p(iV)2 ■ 


(50) 


1 2A2 

c(A^) ^ ap{N) ■ 


(51) 


This solution will be self-consistent with the initial assumption c{N) > 1 
only if ^ 2/0 is negative, which then requires that that /2 and /i should be of 
opposite signs (since Ai is larger than ^)“). Note that for N = 2 the solution 
p(2) exists in the hrst place only if /2 and fi are of opposite signs. Hence 
one concludes that c(2) > 1 is generically expected. Comparing (^) with 
the corresponding equation at the next order iV -|- 1 one obtains at large N 
the recursive relation 


1 

c(A^ + l) 



1 


which can also be written as 

1 1 
c(FTT) ~ c(iV) 



(52) 


(53) 
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showing that 


c(A + 1) < c{N). 

(54) 

Indeed the large N solution of (|5^) is 


1 1 

(55) 

-~ 1 H-. 

c{N) N + K 


with K a constant. Therefore c{N —*• cxo) —»• 1+, that is, though p{N) will 
increase with N, it will approach a constant as iV —> cx). Now, since p{N) 
increases with N, this means that the various approximations that led from 
to will become increasingly accurate. What is remarkable is that 

in the examples studied with c{N) > 1, (^) is already satisfied at low N 
and furthermore the relation too is a reasonable approximation. 

The above relations (|5^|54D were obtained under the assumption c{N) > 
1. If c{N) < 1 then there are apparently no simple relations. For the 
auxiliary series of the beta function in Sec. (3.5), c{N) alternated between 
being slightly larger than one and much smaller than one. 

Consider now the Eq.(^) which is valid at large N and large p{N), 


ASn = 


Ai f 1 
p{N) [c{N) 


( 66 ) 


Remarkably, this simple equation summarizes most of the observed trends. 
When c{N) > 1, (|5B[ ) implies that for Ai < 0, which corresponds to /i < 0 
and minima solutions to (O), ASn > 0, which is indeed observed in the 
examples studied and leads to lower bounds. Similarly if Ai > 0, which 
corresponds to /i > 0 and maxima solutions, ASn < 0 , implying upper 
bounds. The only exception observed so far is the Euler-Heisenberg series 
in Sec.(3.4) which had c{N) > 1, Ai < 0, and yet gave upper bounds. 
Presumably for that example the approximation (|5^) is not appropriate. 
Now suppose that c{N) < 1, as happens roughly for the auxiliary series in 
Sec.(3.5). Then Eq. (|56|) implies that minima, corresponding to /i < 0 and 
hence Ai < 0, give AS' < 0, which explains some of the oddities observed in 


The equation (|5BD also explains the rapid convergence of the bounds. 
Indeed one sees that convergence can be achieved simply by one of two ways. 
Firstly, if c{N) —>■ 1 as N —> oo, that is p{N) —>• po- This type of behaviour, 
is manifested, for example, by the toy-model in Sec.(3.5) of Ref. 0. The 
second way is for c{N) > 1 for all N, which leads to p{N) increasing with N. 
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This is what has been generically observed. Actually, as shown above, the 
second behaviour also gives c{N) —>■ I"*", and hence the convergence is doubly 
fast. Explicitly, from Eqns. (|^ , |5^ , |^) one deduces for the c{N) > 1 case. 




aAx 

aAi 


c(N) 

1 


2 A 2 {N + Kf 


(57) 

( 68 ) 


In summary, for practical problems, the strongest statements can be made 
for the case when c{N) > 1 is observed for the given terms of a partial 
perturbative series, that is when p{N) increases with N. Then the various 
approximations leading to the above equations become increasingly accurate 
at large N, allowing one to make assertions about all N. Firstly, one deduces 
c{N —>■ cx)) —>■ 1’'“. Secondly, for Sn which are global (and local) minima, 
a Sn < Sn +1 for the given terms of the series, the trend will continue for 
larger N, the convergence of the Sn will be rapid, and they will form lower 
bounds to the exact result. However if Sn > Sn+i is observed for the global 
(and local) minima, (as in Sec. (3.4)), then though the trend will continue 
and though the convergence of the Sn will be rapid, it is not a priori obvious 
that they will be upper bounds to the exact result, because then key pieces 
( [T^ - pr5D of the argument are missing. (The same loophole occurs for the 
approximants Sn formed from the p{N)A.) 

If it is observed that c{N) < 1 for the given terms of a series, then the 
various equations above are not necessarily accurate at higher N. In that 
case one can only conjecture that the observed monotonic and rapid conver¬ 
gence will continue at higher orders. 


A2 


Consider now the slope of the approximants Sn{X). From (|4^) , as A 
SnW ^ fo + /lA, so that 


0 , 


dSM 

dX 


A=0 — fl ■ 


(59) 


Thus all the bounds approach the origin with the same value (/o) and slope 
{fl) independent of p and A.This fact can be seen in all the figures. Consider 
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next Eq.(P^ for large p{N), 

dS]\f I 1 Ai 

“ar''” ~ Tap{N) ' 


(60) 


Since the sign of Ai is the same as the sign of /i, this shows that the curves 
S'Ar(A) are not expected to change direction as A varies. This is indeed ob¬ 
served, and explains why the bounds SnW to S'(A) are also good estimates 
of *S'(A) itself only when the latter has a slowly varying Erst derivative. 

As was observed in Sec.(3.5) (see also 0), the approximants 5(A) formed 
from the local extrema p{N) had a more varying slope. In order to understand 
this, set for simplicity /o = 0 so that Sn{0 ) = 0 and let us demand that 


5Ar(Ao) — 0, 


(61) 


SO that SnW curves back to its value at the origin, and so can better ap¬ 
proximate functions like sin(A7r). The condition (|^ then leads from (^3]) 
to 

N A 

0 = V . (62) 

Note that since the p(iV) also have to satisfy the extremum condition (|lT|), 
eqns. (pd] , |62D are actually two coupled equations for Aq and p(N) which we 
would like to analyse for consistency. Firstly, (^) can be used to eliminate 
the leading term in the ’barred’ version of (^), so that now 


dSjv y- (n- l)An 

which at large p{N) gives the slope 

dS]\f . 1 A 2 

“ar''” ~ yp(iv)2 ■ 


(63) 


(64) 


If /2 is opposite in sign to /i then by comparing ( |5^ with ( ^^D one sees that 
indeed the approximants S can change direction Of course this just shows 
that the approximate large p{N) analysis above is self-consistent. Now, in 
Sec.(2), for fi <0, /2 > 0, the p{N) have been defined as positions of local 
maxima when the p{N) are positions of local minima. We can compare the 
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relative magnitudes of the two values as follows. For large p{N), (|H^) has 
the approximate solution 

P{N) ~ . (65) 

By contrast the p{N) are approximate solutions of (|4^) with the left-hand- 
side deleted, which is equivalent to ignoring the mild p dependence of the 

^n’s, 

_o A 

p{N) ~ ^ . (66) 

Thus the values of p{N) are expected to be smaller than those of p{N), and 
the interested reader may verify from the given examples that this is indeed 
the case. Reversing the logic of the argument above, one concludes as follows. 
For that N when p{N) is the position of a local minimum, one expects a local 
maximum at p{N). While the Sn curve is expected to be monotonic in A, 
that of S{X) will not be monotonic if the value of p{N) is smaller than that 
of p{N). 

A similar discussion can be carried out for the curves formed from an 
auxiliary series S'^. Again set for simplicity /o = 0, and dehne 



A 


(67) 


The extremization (1TI|) in p is done with respect to the resummed auxiliary 
series and one deduces for large p (iV), from an equation analogous to 


(|4^), that 


sign 



Ao 


sign(/2). 


( 68 ) 


Now demanding 


*S'Ar(Ao) — 0, 


(69) 


for the reconstructed resummed series S'at = AR^y gives. 


OSn, (dS'^X 


(70) 


which when combined with (| 68 D shows that for fi and /2 of opposite signs, 
the slope of the reconstructed S'at at A = Aq is opposite in sign to its slope 
at the origin which is given by 
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Thus if /i and /2 are of opposite signs, then the auxiliary series may be 
expected to give a reconstructed Sjy with a slope that has variation in sign, 
compared with the slope of the approximant Sn which is obtained by direct 
means, if the two conditions (0 and (|6^ for S have a consistent solution 
p'{N) for some Aq. 
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Figure Captions 

Figure (la): Plots of the divergent series Sp^{\) for the model in Sec.(3.1), 
together with the exact result *S'(A). Starting from the lowest curve and mov¬ 
ing upwards, one has S^, S 3 , S, S 2 , 84 ^. 

Figure (lb): Plots of the resummed series Sn{X) for the model in Sec.(3.1), 
together with the exact result *S'(A). Starting from the lowest curve and mov¬ 
ing upwards, one has S 2 , S 3 , S 4 , S. 

Figure (Ic): Plots of the resummed series 5 '(o) 7 v(A) for the model in 
Sec.(3.1), together with the exact result 5'(A). Starting from the lowest 
curve and moving upwards, one has 5 * 4 , S, S 5 , S 3 . Compared to Fig.(lb), 
the convergence is faster but not monotonic. 

Figure (2a): Plots of the series Sn{X) for the model in Sec.(3.2), together 
with the exact result *S'(A). Starting from the lowest curve and moving up¬ 
wards, one has ^ 3 , S, § 4 , S 2 . 

Figure (2b): Plots of the resummed series Sn{X) for the model in Sec.(3.2), 
together with the exact result 5'(A). Starting from the lowest curve and mov¬ 
ing upwards, one has S 2 , S 3 , S 4 , S. 

Figure (3a): Plots of the resummed perturbation series Sn{\) for the 
model in Sec.(3.3), together with the full perturbative result Spert{X). Start¬ 
ing from the lowest curve and moving upwards, one has S 2 , S 3 , S 4 , S 3 , Spert- 

Figure (3b): Plots of SexactW, S^i^X) and S^i^X) — Snp{X) as dehned in 
Sec.(3.3). At A = 10, the lowest curve is of S^ — Snp, and the highest one is S 3 . 

Figure (4): Plots of the resummed Euler-Heisenberg series Sn{X), to¬ 
gether with Schwinger’s exact result 5'(A). Starting from the lowest curve 
and moving upwards, one has S, Sq, S 4 , S 2 . The approximants approach 
the exact result from above. 

Figure (5a): Plots of the resummed beta function Sn{X) of Sec.(3.5). 
Starting from the lowest curve and moving upwards, one has Ay, Sq, S 3 , S 4 , 83 , 82 . 

Figure (5b): Plots of the resummed series S]\f{X) for beta function in 



Sec.(3.5). Starting from the lowest curve and moving upwards, one has 
S 3 , S's, S'y. The approximants appear to form upper bounds. 

Figure (5c): Plots of the resummed beta function SnW of Sec.(3.5), 
obtained through the auxiliary series. Starting from the lowest curve and 
moving upwards, one has S2, S3, S4, S3, Sq. The approximants appear to 
form upper bounds. The curves for N = 2 and N = 3 are indistinguishable, 
and similarly, those for iV = 4 and N = 5 are very close. 

Figure (5d): Magnihcation of the N = 6 curve of Fig.(3.5) near its non¬ 
trivial zero. 
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